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Initially, we will consider some problems of theoretical mechanics, which are reduced to 
solving ordinary differential equations (ODES). To solve them, the built-in functions of the 
MathCAD package will be used. 


MathCAD function for solving an ordinary differential equation of the m -th order. Recall 
that an ordinary differential equation of the m-th order is called the equation 
y™ = f(ty Y y), (1) 


d™ y(t) 


where y™ = Tt) 


—means an m-th order derivative with respect to the variable t. 


The Koshi problem for such an equation is to find such a function y(t), which is a solution to 
equation (1) and satisfies the initial conditions: 


Hiya? Ce) = Via PU) = Fae (2) 


To solve the Koshi problem considered in the MathCAD package, the odesolve function is 
most often used, which is referred to in the Given decision block. In this block, the m-th order 
differential equation itself and the initial conditions are set. 


The appeal to the odesolve function has the form: 
odesolve(t,t;, n), (3) 
where t - is the name of the variable (the argument of the desired solution); 


tkg — is the value specifying the end of the integration interval (the beginning of the interval 
to); it is specified in the initial conditions set in the Given block by the following 
construction: y = (to) = Yo,y' (to) = yı and etc. (bold characters are used in this case 
<<equals>>); 


n - is an optional parameter that determines the number of intervals into which the segment 
(to, tk) is divided. If this parameter is omitted, then the calculation of n occurs inside the 
odesolve function. 


The result of the work is the value of the solution of the differential equation at time t. 


The use of the odesolve function is shown in the following example. 
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Example 1. The problem of free oscillations of a dynamical system. Let the load Q rest on an 
elastic spring. The deviation of the load from the equilibrium position and at the initial 
moment of time is denoted by yg, and the velocity values are yı (initial conditions). The 
deviation down is positive, up is negative. In the equilibrium position, the weight is balanced 
by the spring elasticity. It is assumed that the force tending to return the load to the 
equilibrium position — the so-called restoring force - is proportional to the deviation, i.e. 
equal to — k'y, where k is some constant value for this spring (<<spring stiffness>>). It is 
proposed that the movement of the load Q is hindered by a resistance force directed in the 
direction opposite to the direction of movement and proportional to the speed of movement of 


the load, i.e. the force -A-v = —A- z, where A = const = 0 (shock absorber). Find the 
function y(t) —the coordinate of the movement of the load. 


Decision. The differential equation of the movement of the load on the spring according to 
Newton 's H law has the form: 


Ë? key. 2 
a= k-y-aA (4) 
where k, A > 0. 


The solution of the Koshi problem is shown in the MathCAD document. 


Qy" (t) = -k-y (t) - A-y' (t) 
y(0)=yoọ y'(0)=0 


Ye `= odesolve(t ,6 , 100) Peres S 


- (0.01 | \ 
Yalt) = € -| 0.01-cos(2-t) +| ——-sin(2-1) | 
\ « J 


Fig.1. Solution of the Koshi problem for example 1. 


All system parameters are set here. To enter a stroke (the designation of the derivative), press 
[Ctrl]+[F7]. The paper presents graphs of the numerical solution (curve y,(t) —solid line) 
and the analytical solution (curve y,(t) —cretics). A good coincidence of the two curves is 
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visible. 


To illustrate the possibilities of MathCAD, consider an example with a spring, the rigidity of 
which varies depending on the sign of the deviation of the load. 


Example 2. The problem of free oscillations of a dynamic system with variable stiffness. In 
the conditions of example 1. we make the following assumption, namely, the stiffness of the 
spring is set as follows: 


_ (100, y(t) >0 
KO) = {ao y(t) <0. ©) 


Decision. To set the variable stiffness (5), we will use the if function. In Fig.2. A document 
solving the Cauchy problem is shown. Comparing the graphs of numerical solutions y,(t) 
(Fig.1.), one can notice the difference between the solutions obtained. The analytical solution 
cannot be obtained, and this graph of the curve y, (t) (Fig. 2.) is missing. 


Yo = 0.01 Ya = 0 


Given 


Q-y"' (t) = -if (y (t) < 0,20,100)-y(t) — A-y' (t) 
y(0)=yo y'(0)=90 


Ye `= Odesolve(t . 6) t := 0.0,0.1..6 


Fig.2. Solution of the Cauchy problem for example 2. 


MathCAD functions for solving systems of first-order differential equations. Many mexanical 
processes they are described by systems of equations containing derivatives of several 
unknown functions for the same independent variable. Consider such a system of equations: 


V1 = filt, Yr Yz, Ym), 
Vo = f(t, Y Yz =» Ym), (6) 


= falt, Ya Yas i nd: 
where y1, V2, ., Ym —are the desired functions of the variable t; 


fa fo, =» fm —set functions of several variables 
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We will consider the Koshi problem for the system (6): 

Y1 (to) = Y1 Y2(to) = Yar» Vm (to) = Ym- (7) 

where y1, V2, -+-» Ym —Set values. 

In this case, the solution is calculated on the interval (to, tg). 

Note that the first-order differential equation is a special case of the system (6) for m = 1. 


The Koshi problem for an ODE system is reduced to solving the Koshi problem for an ODE 
of the form (1). By introducing new variables: 


Yi =Y’ = YY =Y" = Var Im =Y™, (8) 
we come to the first - order ODE system: 
Yi = Y» 
J= Yy 
AENEAN (9) 
Ym-1 = Ym 


Ym = f (6 Yr Var Ym) 
with initial conditions (2): 
To solve first-order ODE systems, the following functions are defined in the package: 
> rkfixed — the Runge-Kutta method with a fixed step; 
> Rkadapt — Runge — Kutta method with adaptive step selection; 
> Bulstoer — the Bulirsch—-Stern method. 
The reference to these functions has the form: 
rkfixed(y,to, tg, n, D); Rkadapt(y, to, tk, n, D); Bulstoer(y, to, tg, n, D). 
The arguments of the functions are the same: 


y - is the vector of initial conditions of dimension m, the projections of which are determined 
by the values (7); 


to, t,-are the boundary points of the interval on which the solution of the differential equation 
is sought; 


n- is the number of points (not counting the initial one) in which the values of the desired 
solution are calculated; 


D- is a vector function of two variables D(x,y) containing a vector of the first derivatives, i.e. 
the right part of the systems (6) or (8). 


The result of these functions is a matrix of size (n + 1) x (n + 1) the first column of which 
is the values of the argument ¢ in which the solution was calculated. The second column and 
the subsequent ones are the numerical values of the solutions y;(t), y2(t), ..., ¥m(t) at points 
to, ty, ...,t, (the first column of the matrix). 


Example 3. Using the rkfixed function to solve the problem Koshi of example 1, reducing it 
to the Koshi problem for a system of first-order ODES. 


Decision. Let's introduce new variables: 
yı (t) = y(t) — the desired solution of the ODE; 
y2(t) = y'(t) —is the first derivative of the solution of the ODE. 
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Then we get a system of first - order ODES: 
yi = Yz 
Pee i (11) 
Y2 5 ~ Gui — gyz 
The initial conditions are y,(t)) = 0,01, y,(t)) = 0. 


The solution of this system is given in the document shown in Figure 3. Note that the 
ORIGIN variable is set to 1 (to match the indexes of variables with the system designations 
(11). The variables y4, yz are elements of the vector, and therefore their index expressions are 
entered after pressing the "left square bracket" key. The variables y,, yz (set the initial 
conditions) are scalar quantities, 


and their indices are entered after pressing the key [-] — decimal point. 


Z := rkfixed (Y „to. tg”, D) 


r =Z” =Z? i=l1.n+1 


Fig. 3. Solution of the Koshi problem for a system of ordinary differential equations 
(example 3) 
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